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We present a numerical approach to simulate the Ferromagnetic Resonance (FMR) of micron and 
nanosized magnetic elements by a micromagnetic finite difference method. In addition to a static 
magnetic field a linearly polarized oscillating magnetic field is utilized to excite and analyze the 
spin wave excitations observed by Ferromagnetic Resonance in the space- and time-domain. Our 
continuous wave approach (CW) provides an alternative to the common simulation method, which 
uses a pulsed excitation of the magnetic system. It directly models conventional FMR-experiments 
and permits the determination of the real and imaginary part of the complex dynamic susceptibility 
without the need of post-processing. Furthermore not only the resonance fields, but also linewidths, 
ellipticity, phase relations and relative intensities of the excited spin wave modes in a spectrum can 
be determined and compared to experimental data. The magnetic responses can be plotted as a 
function of spatial dimensions yielding a detailed visualization of the spin wave modes and their 
localization as a function of external magnetic field and frequency. This is illustrated for the case 
of a magnetic micron sized stripe. 


I. INTRODUCTION 

The detailed understanding of spin wave spectra of 
magnetic micro- and nanostructures and their magneti¬ 
zation dynamics has found increasing interest from both 
fundamental and applied points of view for example in 
spin caloritronics and spin torque phenomena 1 4 . A pow¬ 
erful tool to investigate these spin wave spectra exper¬ 
imentally is the Ferromagnetic Resonance (FMR) de¬ 
tected in the frequency domain 5,6 . However in most cases 
the obtained FMR-spectra are complex in nature fea¬ 
turing several -often overlapping- resonances and require 
theoretical descriptions of the nanostructured magnetic 
systems to extract quantitative information. Micromag¬ 
netic simulations of the FMR can be used to model those 
systems and provide additional information on the char¬ 
acter of the observed magnetic excitations as well as their 
dependence on magnetic parameters, geometries, confine¬ 
ment effects or charge currentss 7,8 . This is especially of 
interest when the complex geometries and interactions of 
the nanoscale ferromagnet aggravate quantitative analyt¬ 
ical approaches. 

Here we present a finite difference method utilizing a 
homogeneous oscillating magnetic field to simulate FMR 
spectra corresponding to experiments. In addition we 
show how to further analyze the spectra by visualizing 
the spatial distribution of the magnetic excitations. We 
start by describing the problem definition, initialization 
and recorded data of the simulations in section II. Sub¬ 
sequently in section III the derivation of FMR-spectra is 
described in detail as well as determination of the FMR 
fields, linewidth and ellipticity of the resonances. In sec¬ 


tion IV we investigate the resonances contained in the 
spectra in terms of spin waves and spatial variations. 

II. METHOD 

The micromagnetic simulations presented here, are 
based on the public domain 3D-OOMMF (Object Ori¬ 
ented MicroMagnetic Framework ) 9 solver. This finite dif¬ 
ference software solves numerically the Landau-Liftshitz 
equation 15 (LLE) 

“f = -y (MX ft.) -2^ (MX (MX A,)) 

( 1 ) 

where 7 is the gyromagnetic ratio, a the Gilbert damp¬ 
ing constant and H e ff the effective field. As time evolver 
a Runge-Kutta method is used. Further details on the 
implement ion can be found in Ref. 5. Our approach to 
simulate FMR-experiments can be split into three differ¬ 
ent steps: 1. relaxation 2. transient phase 3. dynamic 
equilibrium. 

For initialization the spatial dimensions of a nanos¬ 
tructured ferromagnetic system are defined by a grid of 
equal rectangular cells. All cells are assigned with an 
identical magnetization vector M , located at the cen¬ 
ter of the cell. In order to simulate the external field of 
FMR-experiments, a static magnetic field H is applied to 
the system. To obtain the static magnetic ground state 
dM /dt = 0 a relaxation simulation is performed, without 
applying any excitation. So the motion of M damps out 
and M will reorient to an equilibrium direction, given 
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by the local effective field. For faster convergence the 
precession term in the LLE may be switched off and the 
damping constant a set large. As stopping criterion for 
the simulation typically values of du/dt < 0.001 °/ns (u 
is the unit vector of M ) are chosen, to achieve the quasi 
static state, which is used as the inital state for the sub¬ 
sequent FMR-simulations. 

In addition to the static field H a linearly polar¬ 
ized oscillating magnetic field (t) = ho • sin (uj • t) 
is added for continous excitation of the magnetization 
(CW). Hrf is uniform over all cells, oriented perpendic¬ 
ular to the static field and corresponds to the magnetic 
microwave field used in conventional FMR-experiments: 
ho = 398 A/m = 0.5 mT (satifying the relation ho « 
H). 

Due to the interaction with Hrf a driving torque is ex¬ 
erted on the magnetization. After a transient phase the 
magnetization reaches a dynamic equilibrium, precess- 
ing around the effective field with the angular frequency 
uj. In this state Krf transfers power to the magnetic 
system to compensate dissipation, induced by damping. 
To study the dynamic equilibrium and discard transient 
effects, a fixed time period is simulated, without gen¬ 
erating data for analysis. This time period is given by 
t s =27 r • N/uj , using an integer oscillation number N of 
Jirf (details described in section III). N is typically set 
between 40 and 60 depending on the magnetic parame¬ 
ters (e.g. damping constant a). This provides a constant 
precession amplitude between consecutive oscillation cy¬ 
cles of the magnetization with a deviation of less than 
0.02 %. When the simulation time reaches t s the actual 
parameters (like magnetization M , oscillating field Jirf , 
static field H ,...) are stored for each following iteration 
step of the time evolver. This process continues for one 
further cycle of the oscillating field (N + 1) and consists 
of at least 1000 iteration steps, which provides a time 
resolution in the ps regime. By analyzing these data, it 
is possible to follow the precession trajectory of the cell 
specific magnetization in the space- and time-domain. 

III. ANALYSIS OF THE SIMULATION 

We now describe the analysis of the simulated data for 
the dynamic equilibrium. In the linear response regime 
discussed here the dynamic magnetization m (t) of each 
cell exposed to the oscillating field is described by the 
dynamic susceptibility-tensor % = y' + i y": 

fn(t) = x hRF (t) = X h 0 exp -i(wt - n/2) (2) 

Where we have used a complex representation of the 
susceptibility and the applied sinusoidal varying mag¬ 
netic field Jirf • % is a 3 x 3-tensor with elements x ij • 

To illustrate the motion of m in respect to the hRF , 
the simulated time dependences for the case of an infinite 



FIG. 1. Micromagnetic simulation of the time dependent re¬ 
sponse of the dynamic magnetization driven by an external 
dynamic field hRF for an infinite thin film in the xy-plane. 
The normalized dynamic magnetization in the yz-plane (solid 
lines) are shown together with Hrf (dashed line) oriented 
along the film-normal in the z-direction. The static field H is 
oriented in the film plane (x-direction) and is chosen to match 
the resonance condition. 

film spanning the xy-plane are shown in fig. 1 . Jirf is 
oriented in the z-direction (out of plane) and H in the 
x-direction (in plane), respectively. The displayed time 
dependent components of the oscillating dynamic mag¬ 
netization (solid lines) lie in the yz-plane driven by hRF 
(dashed line). The oscillation of m and its components 
rriy iZ is described by its amplitude frequency uj and 
phase relation <j) y , z in respect to hRF . Note that the 
driven component m z exhibits a phase-shift of 90° to the 
oscillating field as expected for a resonantly driven sys¬ 
tem. The precessional motion of the magnetization in 
equilibrium as well as the ellipticity of its trajectory can 
be readily observed by the 90° phase shift between the 
dynamic components m yiZ and the ratio of their differing 
maximal amplitudes. 

The simulated FMR-spectra, which can be quantita¬ 
tively compared to experimental ones are derived as de¬ 
scribed in the following. The power P absorbed by the 
magnetic system from hRF and consequently the FMR- 
signal S is proportional to the diagonal element x zz °f 
the imaginary part of the susceptibility 16 : 

s OC P oc Xzz (3) 

To determine x'zz f rom the simulation it is sufficient 
to analyze the component of m parallel to hRF (in this 
case m z ) after a complete cycle of oscillations as given for 
the time t s (see also section 2). This can be seen when 
considering the observable real part of fh in equation 2: 

m z = 5R( Xzz h 0 exp(-i(ut-n/2))) (4) 

Inserting t s (u)t s = 2 tt N) yields: 
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m z (t s )= - x'L h o (5) 

Hence, a proportional FMR-Signal (x zz ) can s ^ m_ 
ulated by directly monitoring — m z (t s )/ ho without the 
need for extensive post-processing. (By a similar logic) 
the real part of the susceptibility x'zz ~ an d therefore the 
complete Xzz - can be retrieved from the simulation by 
extracting m z at the maximum of Hrf • 

Fig. 2 shows the simulated external field dependent 
amplitude, phase and imaginary part of the susceptibil¬ 
ity for an infinite thin film with the magnetic parameters: 
exchange constant A = 13 • 10 -12 J/m 3 , saturation mag¬ 
netization M s = 8.3 • 10 6 A/m, g-factor g = 2 . 12 , Gilbert 
damping parameter a = 5 • 10 -3 . 

The magnetic response shows the typical hallmarks of 
a driven oscillator in respect to phase and amplitude and 
a lorentzian absorption curve 17 . This enables one to de¬ 
termine the resonance positions as well as their linewidth 
and relative signal strength. In the case of multiple res¬ 
onances a decomposition of the resultant spectra (exper¬ 
iment or simulation) into lorentzian absorption lines is 
needed to obtain those quantities. Subsequently this can 
be compared to complex experimental results as for ex¬ 
ample for the case of magnetic micron sized stripes 11 ^ 13 . 



FIG. 2. Micromagnetic simulation of amplitude, phase and 
imaginary part of the susceptibility for an infinite thin film 
for different static fields and a driving frequency of 15 GHz. 
The magnetic resonance exhibits the classical hallmarks of a 
driven oscillator, where the imaginary part is proportional to 
the experimentally detected FMR-Signal. 


IV. CHARACTERIZING THE MAGNETIC 
EXCITATIONS 

To further investigate the magnetic resonances in the 
simulated FMR spectra the amplitude, phase and imag¬ 
inary part of the susceptibility may be analyzed for each 
cell of the magnetic system individually. The retrieved 
spatial distribution of the magnetic response often helps 


to identify spinwave like resonances, their wavevector or a 
localized character of the excitations. Here we exemplar- 
ily simulate a FMR Spectrum of a 5 pm x 1 pm x 20 nm 
stripe at 14 GHz. As magnetic parameters we choose: 
exchange constant A = 13 • 10 -12 J/m 3 , saturation mag¬ 
netization M s = 8.3 • 10 6 A/m, g-factor g = 2.13, Gilbert 
damping parameter am 7 • 10 -3 . The simulated FMR- 
Signal together with the stripe geometry and field orien¬ 
tation is shown in fig. 3. 



FIG. 3. Micromagnetic simulation of the susceptibility x'zz 
for a 5 pm x 1 pm x 20 nm Permalloy stripe and a frequency of 
14 GHz. In the confined system multiple magnetic resonances 
are observed (labeled from 1 to 4) 

In such a confined system multiple magnetic reso¬ 
nances (labeled from 1 to 4) with differing positions, 
linewidths and intensities occur 18 . The spatial distribu¬ 
tion of the normalized imaginary part of the susceptibil¬ 
ity for the most intense resonance 1 is shown color coded 
in fig. 4. As can be seen this magnetic resonance ex¬ 
hibits the strongest excitation in the center of the stripe 
and will here be referred to as a localized quasi-uniform 
mode. A different mode character can be observed for the 
less intense magnetic resonance 2 in fig. 4. x" zz shows 
a change in sign across the stripe, two nodal lines and a 
wavelike varying dependence along the dashed line. 

This is very well approximated by a sinusoidal function 
as shown in fig. 5 and resembles the expected modepro- 
file of a standing spinwave with wavelength 734 nm and 
dipolar pinning 19 at the edges of the stripe, which devi¬ 
ates from simply closed or open pinning conditions. In 
a similar consideration the resonances 3, 4 can be as¬ 
signed to higher standing spinwaves across the width of 
the stripe, where the wavelength decreases for smaller 
resonance fields. 

By such a spatial analysis of the magnetic excitations 
one can for example explore the dependence of the reso¬ 
nance position, linewidth, mode profile, ellipticity, inten¬ 
sity and pinning conditions on the magnetic parameters 
as well as on the exact geometry of the magnetic systems. 
This information can be crucial for planning experiments 
and lead to a deeper understanding of the multiple reso- 
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FIG. 4. Normalized imaginary part of the susceptibility for 
the magnetic resonances labeled as 1 and 2 respectively. Res¬ 
onance 1 shows a localized character (quasi uniform mode) 
with the strongest excitation in the center region of the stripe. 
Resonance 2 exhibits a harmonic dependence across the width 
of the stripe as expected for a standing spinwave with a wave¬ 
length of 734 nm. The profile along the dashed line is shown 
in fig. 5. 


nances in experimentally observed spectra. 



y [|jm] for x = 2.5 |jm 

FIG. 5. Lineprofile of the normalized imaginary part of the 
susceptibility (black dots) for the magnetic resonances labeled 
as 2 along the dashed line of fig. 4. The spatial variation 
is very well approximated by a harmonic dependence of a 
standing spinwave with a wavelength of 734nm (red line). 


V. CONCLUSION 


In summary, we presented a numerical approach to 
simulate the FMR of nanosized elements and nanostruc- 
tured systems using a finite difference method. The sim¬ 
ulation utilizes a homogeneous linearly polarized oscil¬ 
lating magnetic field to drive the magnetic system. We 
would like to point out that this method is very easy to 
implement and handle. The derivation of simulated FMR 
spectra, which directly correspond to experimentally de¬ 
tected spectra, is explained in detail. This allows one 
to compare not only the resonance position, but also the 
linewidth and intensities of experiment and simulation 
directly. 


Furthermore the spatial dependence of the magnetic 
excitations spin waves and localized resonances can be 
identified and visualized. Additionally the time depen¬ 
dent trajectory of the magnetization as well as the accu¬ 
rate phase relation to the driving field can be extracted. 
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